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Abstract 

We present here some newer characteristics pertaining to paramagnetic Meissner effect like response in a single crystal of the 
low Tc superconducting compound 2H-NbSe2 via a detailed study of effects of perturbation on the field-cooled magnetization 
response. In the temperature range, where an anomalous paramagnetic magnetization occurs, the field-cooled magnetization 
response is found to be highly metastable: it displays a curious tendency to switch randomly from a given paramagnetic value 
to a diamagnetic or to a different paramagnetic value, when the system is perturbed by an impulse of an externally applied ac 
field. The new facets revealed in a single crystal of 2H-NbSe2 surprisingly bear a marked resemblance with the characteristics of 
magnetization behaviour anticipated for the giant vortex states with multiple flux quanta (I/4>o, 'f’o = hc/2e, L > 1) predicted 
to occur in mesoscopic-sized superconducting specimen and possible transitions amongst such states. 

PACS numbers: 74.25.Ha, 74.25.Op 
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I. INTRODUCTION 


The observation of an anomalous paramagnetic mag¬ 
netization signal (viz., paramagnetic Meissner effect 
(PME)), instead of the usual diamagnetic behaviour, on 
field - cooliM a superconducting specimen, continues to 
attract [l|, Q attention ever since its discovery Q in ce¬ 
ramic sample(s ) of a high Tc cuprate. In addition to high 
Tc cuprates pMl^. the PME-like attribute is known to 
occur in a wide variety of other superconductors [31^. 
Numerous explanations (see Ref. [27| for a review), e.g ., 
the d - wave superconductivity @ , orbital glass [13, , 

presence of tt contacts leading to spontaneous currents 
H, H, Josephson junctions [s^, etc., have been ad¬ 
vanced as the origin of PME like signal(s) in high Tc 
superconductors. A model proposed by Koshelev and 
Larkin [3l| envisaged the possibility of trapping of mag¬ 
netic flux in the interior of a superconducting sample due 
to inhomogeneous cooling. In such a circumstance, the 
PME can occur due to unbalancing between two oppo¬ 
sitely directed currents, viz., (i) a (paramagnetic) cur¬ 
rent set up in the interior that attempts to screen the 
effect of trapped flux and, (ii) a usual (diamagnetic) cur¬ 
rent flowing on the surface that shields the external field 
[ 31 II . However, within the Ginzburg - Landau framework, 
the theoretical works by Moshchalkov et al. [s^l and 
Zharkov [sl], for model cases of mesoscopic cylindrical 
shaped superconductors, lead to the proposition that the 
multi-quanta vortex matter, i.e., vortices with multiple 
flux quanta, L^o, <i)o = hc/2e, L > 1, nucleating be- 
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low a third critical held (Hcs), at the onset of surface 
superconductivity can also give rise to the PME. 

Recently, a controlled switching of the PME into the 
usual (diamagnetic) Meissner effect has been vividly 
demonstrated by Xing et al. [ 2 ^ in a superconducting 
ferromagnet Pb-Go nanocomposite, wherein the source 
of PME is argued to be related to some different mecha¬ 
nism other than those stated above. This manipulation 
of PME is possible by the change in orientation of mag¬ 
netic moments of Co nanoparticles with respect to the 
external magnetic field [ 2 ^ . Such a tuning of the anoma¬ 
lous PME and the non-anomalous (diamagnetic) Meiss¬ 
ner effect has rarely been reported in the literature. As 
another case of manipulating a PME like response, we 
have explored some novel facets of the PME in a single 
crystal of 2H-NbSe2 [s^, an example of a conventional 
low Tc superconducting system often studied to explore 
novel notions in vortex matter. The results being re¬ 
ported ahead in 2 H-NbSe 2 present an advancement over 
studies earlier reported by some of us in Ca 3 Rh 4 Sni 3 
and Nb [i3- The curious behaviour pertains to an un¬ 
usual paramagnetic magnetization (a la PME) and its 
manipulation just below Tc in the (isofield) temperature- 
dependent dc magnetization {M{T)) data. We And that 
in a short window of temperature, where the PME oc¬ 
curs below Tc in 2H-NbSe2, an external perturbation can 
randomly change the magnetization from a given param¬ 
agnetic value to a larger or lesser paramagnetic value or 
from a given paramagnetic state to a diamagnetic one. 
Across the same temperature interval, the system can 
also transit from a given diamagnetic state into a param¬ 
agnetic one after being perturbed by an ac held impulse. 
The vortex matter in the field-temperature domain of the 
PME in 2H-NbSe2 is thus found to be highly metastable. 
The results in a bulk 2H-NbSe2 specimen unexpectedly 
appear to echo the consequences of nucleation of several 


1 



2 


III RESULTS 




Figure 1: (Color online) (a) Temperature variation of Mzfc, 
Mfcc and Mfcw in a field oi H (|| ab) = 50 Oe. The inset in 
panel (a) illustrates Mzfc (T) plot in the range, 1.8K < T < 
Tc. Panel (b) elucidates M{T) plots on the expanded scale 
in the vicinity of Tc. Anomalous paramagnetic signals are 
evident in all three modes. The paramagnetic magnetization 
reaches its extreme limit at Tp during the field-cooled cool¬ 
down run, Mfcc (T). 

metastable multiple flux quanta states (L > 1) prop osed 
in the context of mesoscopic samples [H, HI, Isfil - lioj l . 


II. EXPERIMENTAL DETAILS 

The single crystal of hexagonal 2H-NbSe2 chosen for 
the present study is platelet-shaped with a planar area 
of 4.12mm^, thickness of 0.17mm and mass of 4.64mg. 
The superconducting transition temperature (Tc) for this 
crystal is measured to be about 7.15K. The crystallo¬ 
graphic ab plane of the hexagonal 2H-NbSe2 coincides 
with the plane of the platelet, while the direction nor¬ 
mal to it corresponds to the crystallographic c - axis 
of 2H-NbSe2. The magnetization measurements were 
performed using a Superconducting Quantum Interfer¬ 
ence Device-Vibrating Sample Magnetometer (SQUID- 
VSM, Quantum - Design, Inc., USA) in such a way that 
the magnetic field was directed nominally parallel to the 
plane of the platelet (i.e., H\\ab), with less than 5 de¬ 


gree error in the angle of alignment. An advantage of 
investigating the sample along this orientation is that 
the demagnetization factor is small ('^ 10“^) and the as¬ 
sociated geometrical/boundary effects are expected to be 
minimal. We estimated, using the flux gate option of the 
SQUID-VSM, the field-inhomogeneity of the supercon¬ 
ducting magnet over a scan length of 8 mm. It is found to 
be only of the order 10“^ Oe. During the magnetization 
measurements, the amplitude of vibration of the sample 
was kept small (« 0.5 mm) in most of the runs (unless 
sgecifically stated otherwise) so that the possible artefact 
due to field inhomogeneity along the scan length of 
the superconducting sample in SQUID - VSM is so minus¬ 
cule that it is of no significance. The effect of changing 
the amplitude of sample vibration on the magnetization 
value (at a given H, T value) in the superconducting state 
of a standard Indium sample has also been investigated. 
We registered no difference in the magnetization response 
of Indium specimen in its superconducting state, when 
the amplitude of vibration in VSM measurements was 
changed from 0.5 to 8 mm. Hence, under the normal 
circumstances, it can be stated that the change in am¬ 
plitude of sample vibration (upto a maximum of 8 mm 
in the present VSM instrument) is not expected to alter 
the dc magnetization value of a superconducting sample 
at a given H and T. In the case of 2H-NbSe2 crystal, we 
have observed a wide variation in (isofield) magnetization 
response, ranging from paramagnetic to diamagnetic val¬ 
ues at different amplitudes (see Fig. 7), which is argued 
to be a consequence of metastability effects prevailing in 
this system, a key focus of the present report. 


III. RESULTS 

A. Paramagnetic signal on field-cooling 

Figure 1 displays the isofield magnetization responses 
obtained in the zero-field cooled {Mzfc {T)), field- 
cooled cool-down {Mfcc {T)) and field-cooled warm-up 
{Mfcw {T)) modes in 2H-NbSe2 in a field of 50 Oe. For 
Mzfc {T) run, the crystal was initially cooled down to 
1.8 K in nominal zero field. A field, H = 50 Oe, was 
then applied (such that H || ab) and the magnetization 
obtained while warming to higher temperatures, yielding 
Mzfc{T). Without changing the magnetic field, the 
M{T) data were again obtained while cooling the sam¬ 
ple back to the lowest temperature {Mfcc {T)). There¬ 
after, the magnetization values were measured again 
while warming it to yield Mfcw {T). The Mzfc{T) 
curve shows the usual fall in diamagnetic response (cf. 
inset panel of Fig. 1(a)) of a superconducting specimen 
as the temperature is swept up. 

An expanded portion of M{T) plots of Fig. 1(a), in 
the proximity of Tc is displayed in Fig. 1(b). Unusual 
paramagnetic responses below Tc can be observed in all 
the three modes. It is pertinent to note that the param¬ 
agnetic signals in Mzfc {T) and Mfcw {T) runs are re- 
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stricted to a very narrow temperature interval (~ 50 mK) 
a little below T^. During Mpcc (T) run, the paramag¬ 
netic response nucleating at Tc enhances as the temper¬ 
ature is lowered towards a temperature marked as Tp. 
Thereafter, the Mpcc {T) values start to decrease, and 
cross over to the diamagnetic values at a temperature of 
about 6.75 K. The diamagnetic Mpcc (T) signal reaches 
a saturated value below 4K (cf. main panel of Fig. 1(a)). 
During warm-up run, saturated Mpcw (T) values re¬ 
main sustained upto about 6K. Above this tempera¬ 
ture, the diamagnetic Mpcw (T) values decrease, and 
they crossover to the paramagnetic values to merge into 
the Mzpc (T) curve at about 7.08 K, which is nearly the 
same temperature (marked as Tp) below which paramag¬ 
netic signal starts to decrease during the corresponding 
Mpcc (T) run. It can be noted in Fig. 1(b), that the 
paramagnetic response in Mpcc (T) run is much larger 
than that in Mzpc (T) run or Mpcw (T) run. This ob¬ 
servation implies that the positive magnetization signal 
below Tc is multi-valued (and thus metastable) and de¬ 
pends significantly on the thermomagnetic history of the 
sample. 

Figure 2 illustrates the Mpcc (T) curves for H (|| ab) = 
100 Oe, 500 Oe and IkOe, which intersect unexpectedly 
at a temperature, marked as Tcross (~ 6.84 K), where the 
magnetization incidentally also crosses the M = 0 axis. 
Below Tcross 1 the magnetic behaviour displayed in Fig. 2 
is such that the diamagnetic response decreases as the 
field increases {{AM/AH)> 0), which is typical of a 
vortex state of a type - II superconductor in the field in¬ 
terval, Hci < H < Hc 2 - On the other hand, the mag¬ 
netization response curves above Tcross imply that the 
paramagnetism decreases as the applied field increases 
|32]. 



T(K) 


Figure 2: (Color online) Portions of the Mpcc (T) curves at 
H — 100 Oe, 500 Oe and 1 kOe, intersecting unusually at a 
characteristic temperature Tcross- 



Figure 3: (Color online) Mpcc (T) data recorded a.t H — 
50 Oe in several runs. In each of these runs, an impulse of 
an ac field (peak amplitude = 10 Oe, frequency = 211 Hz) was 
applied a.t T — 7.07K (encircled), after which the field-cooled 
magnetization change randomly to either paramagnetic or 
diamagnetic values, and Mpcc (T) traverses different paths 
as displayed for three different runs, viz., perturbed (#1), 
perturbed (#2) and perturbed (#3). The pristine Mpcc (T) 
(open circles) was recorded without perturbing the system 
throughout the cooling. The inset displays Mpcc (T) for 
pristine and perturbed (#1, #2, and ffS) cases in the full 
range, 2 K < T < 7.5 K. 


B. Metastability and non - uniqueness in magne¬ 
tization response 

The multi-valuedness in M{T) observed in Fig. 1(b) 
motivated further exploration of magnetization response 
in 2II-NbSe2. We illustrate in the main panel of 
Fig. 3, portions of Mpcc {T) curves in the vicinity of 
Tc recorded at 77 = 50 Oe during several runs where a 
perturbation had been induced at an intermediate stage. 
First, we recorded the magnetization data (see open cir¬ 
cles in Fig. 3) while field-cooling the sample to 2 K, sim¬ 
ilar to the way Mpcc {T) data of Fig. 1 were obtained. 
This curve will henceforth be called as “pristine” in the 
discussion ahead. We then recorded the Mpcc {T) data 
(without changing the external field) while cooling the 
sample down to (near) peak temperature {Tp) of the 
PME, i.e., T Ri 7.07 K, which is encircled in the main 
panel of Fig. 3. At this temperature, the Mpcc{T) 
measurements were paused and the system was momen¬ 
tarily subjected to a perturbation in the form of an ac 
field impulse of peak amplitude 10 Oe (frequency of ac 
field = 211 Hz) imposed for a duration of 6 seconds. It 
was checked that the results being reported in Fig. 3 
do not depend on the duration of the impulse. The 
Mpcc {T) data recorded after the impulse treatment 
showed an abrupt switching from paramagnetic to dia- 
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magnetic value, which is indicated by an arrow point¬ 
ing towards open triangles in main panel of Fig. 3. The 
Mpcc (T) curve (open triangles labeled as “perturbed 
(#1)”) thereafter traverses a path which is akin to the 
one that seems to display the usual (diamagnetic) Meiss¬ 
ner effect in contrast to the anomalous PME observed 
for the pristine curve (open circles). In the next run, 
we again cooled the sample down to the same temper¬ 
ature, Tp (~ 7.07K), and once again applied the same 
impulse. Surprisingly, the Mpcc (T) data now displayed 
(see perturbed (#2) in Fig. 3) an enhancement in the 
paramagnetic signal (shown by upward arrow pointing 
to open squares in Fig. 3). This is in sharp contrast 
to the impact of an impulse in the “perturbed (#1)” 
case, which involved switching of the magnetization from 
paramagnetic to a diamagnetic value. On further lower¬ 
ing of temperature, the Mpcc {T) curve for perturbed 
(#2) can be seen to crossover to diamagnetic values at 
a temperature (« 6.48 K), which is lower than the corre¬ 
sponding crossover temperature of the pristine Mpcc (T) 
curve. In the identical imposition of an ac impulse in 
the third run, the Mpcc {T) data (perturbed (#3)) can 
be switched from paramagnetic to a diamagnetic value 
which is even larger than that during perturbed (#1). 
On progressively lowering the temperature below Tp, the 
Mpcc (T) curve for “perturbed (#3)” (closed squares) is 
seen to stay more diamagnetic (at a given T) than that 
observed in the “perturbed (#1)”, as shown in Fig. 3. 

In the inset panel of Fig. 3, we show the Mpcc {T) 
curves for the pristine and perturbed states (#1, #2 and 
^3) in the entire temperature range of investigation, i.e., 
2K< T < 7.5K. It is curious to note that the satu¬ 
rated 47rM value at 2 K in a given perturbed run roughly 
amounts to the summation of the saturated 47rM for the 
pristine curve and the change in magnetization, A (47rM) 
induced at Tp due to an impulse treatment. 

Subsequent efforts of applying the impulse not only at 
Tp but also in its vicinity, viz., across 6.75 K< T < Tc 
(all data not shown here), lead to the inference that 
the Mpcc (T) value, after perturbation in the PME re¬ 
gion, can unpredictably change very widely from para¬ 
magnetic to diamagnetic or from a given paramagnetic 
to a less/more paramagnetic value. We also applied im¬ 
pulse with different amplitudes and at different frequen¬ 
cies of the ac field, however, the results turned out to be 
same. In all, we can surmise that the vortex matter in 
the temperature regime of PME in 2H-NbSe2 is highly 
metastable and the magnetization associated with such 
metastable states can show a rich diversity, as is evi¬ 
denced in Fig. 3. When subjected to a perturbation, the 
system can transit between these metastable states, as is 
reflected from abrupt change(s) in magnetization. 


C. Path dependent M{T) response 

We now demonstrate in Fig. 4, the impact of an ac field 
impulse applied at various temperatures, during a given 
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Figure 4: (Color online) Mpcc (T) recorded at H = 50 Oe 
in various runs perturbed at different temperatures, (a) Per¬ 
turbed (#1) shows an enhanced paramagnetism at 7K and 
an enhanced diamagnetism at 6K when compared with the 
pristine Mpcc (T). In case of perturbed (#2), an impulse at 
6 K displays equal change in magnetization as is seen in case 
of perturbed (#1). Inset shows pristine and perturbed (#1) 
in the full temperature range, (b) Mpcc (T) traced while an 
impulse is applied at short temperature interval in two runs, 
perturbed (#1) and perturbed (#2), leading to a character¬ 
istic path dependence above 6.75 K, whereas below it, both 
follows the same path. 


Mpcc (T) run at = 50 Oe. The pristine Mpcc (T) 
shown in Fig. 4 is the same as presented in Fig. 3. 
When an impulse was applied at T = 7K (encircled in 
Fig. 4(a)), an enhanced paramagnetic signal was seen 
as depicted by the Mpcc (T) curve named “perturbed 
(#1)”. This curve crosses over to the diamagnetic re¬ 
gion smoothly at T ~ 6.4 K, and thereafter remains less 
diamagnetic (than the pristine Mpcc (T)) until about 
T = 6K, echoing the “perturbed (#2)” case of Fig. 3. At 
T = 6K, an impulse was imposed again (see the upper 
encircling done corresponding to 6 K in Fig. 4(a)) and the 
Mpcc curve thereafter (i.e., perturbed (#1)) becomes 
more diamagnetic than the pristine one. The field-cooled 
curve (perturbed (#1)) follows a path shown by a contin¬ 
uous line below 6 K in Fig. 4(a),which is different from its 
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pristine counterpart. The saturated Mpcc (T) value at 
2 K corresponding to perturbed (#1) curve is found to be 
nearly twice the saturated Mpcc (T) in the pristine case, 
as is evident in the inset panel of Fig. 4(a). In another 
run, we applied an impulse once again at 6 K, but on the 
pristine field-cooled (FC) state (see another encircling 
done at 6K in Fig. 4(a)). Such an impulse treatment 
resulted in the diamagnetic value of pristine Mpcc (T) 
to switch over to a still larger diamagnetic value. The 
Mpcc (T) curve thereafter, i.e., perturbed (#2) shown 
by the dashed line in Fig. 4(a), can be seen to be even 
more diamagnetic than the other two curves. We also 
found that after several iterations, that the change in 
magnetization induced due to impulse treatments below 
about 6.7K is only of one kind, i.e., from a given dia¬ 
magnetic value to another diamagnetic value, unlike the 
behaviour noticed above 6.75K (i.e., in the PME region), 
where one could see switching from para/diamagnetic to 
more or less para/diamagnetic value (see Fig. 3). 

In Fig. 4(b), we show the Mpcc (T) curves at H = 
50 Oe, when an ac field impulse is consecutively im¬ 
posed at short temperature interval(s) in a given FC 
run. We first applied an impulse at T = 7K (encir¬ 
cled in Fig. 4(b)), and observed a switching of param¬ 
agnetic Mpcc of the pristine state into a diamagnetic 
value (perturbed (#1)) as shown by continuous line. 
This Mpcc (T) curve (perturbed (#1)) traces the record 
while cooling the sample further down from 7 K to 6.95 K, 
where the application of another impulse switches the 
diamagnetic value into a paramagnetic value, which is 
even more paramagnetic than the pristine one. As the 
cooling was progressed to T = 6.9 K, an impulse treat¬ 
ment there lead to the change from paramagnetic Mpcc 
value to near zero magnetization value. This process of 
applying an ac impulse at short intervals was continued, 
and the Mpcc {T) curve (#1) displayed in Fig. 4(b) was 
traced. In another field-cooling run, we applied the im¬ 
pulse for the first time at T = 6.9 K, which resulted in 
switching from paramagnetic to a diamagnetic value as 
shown by the dashed curve (perturbed (#2)) in Fig. 4(b). 
This curve was recorded while further cooling, it also 
comprised impulse treatment (s) at several intermediate 
temperatures (see steps in magnetization where the im¬ 
pulse was applied) during the cool-down process. 

An examination of Mpcc (T) data recorded in the 
two runs, “perturbed (#1)” and “perturbed (#2)”, 
in Fig. 4(b) shows a pronounced path-dependence in 
Mpcc (T) at temperatures above about 6.7 K. This sug¬ 
gests that the impact of an impulse above 6.7 K is com¬ 
pletely unpredictable: the system can display any magne¬ 
tization value ranging from paramagnetic to diamagnetic. 
The impulse seems to be acting as a drive which brings 
out transition from a given (paramagnetic/diamagnetic) 
state of the system into another and hence leads to a 
non-uniqueness in Mpcc (T). On the other hand, the 
two Mpcc (T) runs (“perturbed (#1)” and “perturbed 
(#2)”) tend to overlap when an impulse was applied be¬ 
low about 6.7 K. Below this temperature, the two curves 



Figure 5: (Color online) Mpcc (T) in pristine and perturbed 
cases at H = 50Oe. In perturbed (#1), an impulse was first 
imposed at T = 6.7K (encircled), followed by another at 6K, 
and at every IK difference below 6K. Saturated Mpcc (T) 
value of pristine and perturbed (#1) are much different. In 
contrast, the impulse applied at 4K (perturbed (#2)) and at 
3K (perturbed (#3)) brings out minimal changes. No change 
is registered when the impulse is applied at 2 K. 

respond identically to any further impulse treatment and 
the magnetization values remain diamagnetic. 

To explore the highest saturated value of diamagnetic 
response, we chose to successively apply an ac field im¬ 
pulse below 6.7 K at every 1K interval beginning from the 
pristine field-cooled (FC) state created at T « 6.7K as 
shown in Fig. 5. The arrows in Fig. 5 indicate the tem¬ 
peratures at which the impulse was applied. The first 
impulse applied at 6.7 K changes the magnetization from 
a diamagnetic value (encircled) of pristine Mpcc (T) to 
a larger diamagnetic value (perturbed (#1)). After low¬ 
ering the temperature to 6 K, another impulse drives the 
magnetization into a still more diamagnetic state and so 
on. The saturated 47rM value for the perturbed (#1) 
Mpcc (T) curve in Fig. 5 at 2K happens to be about 
— lOOe, which is more than twice the saturated 47rM 
value (ss —4 Oe) of the pristine curve. 

The commencement of application of an impulse at 
temperatures lower than 4 K resulted in a small change in 
the magnetization. For example, the changes in magne¬ 
tization after the impulse treatment at T = 4 K and 3 K 
in pristine Mpcc (T) runs were comparatively smaller in 
the perturbed states “#2” and “#3”, respectively. Even¬ 
tually, at T = 2 K (perturbed (#4)), no change in pristine 
Mpcc was witnessed after the application of an impulse 
at this temperature (cf. Fig. 5). 

The results in Figs. 4 and 5 focus on Mpcc {T) data, 
while subjecting the sample to perturbations at interme¬ 
diate temperatures during cool-down cycle. In Fig. 6, 
we show another set of magnetization data, wherein we 
present the Mpcw (T) responses recorded while sweep¬ 
ing the temperature up after an impulse treatment to 
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(pristine) FC states created at various temperatures. In 
the first run, after getting an enhanced paramagnetism 
for impulse applied at T = 7 K (encircled in main panel of 
Fig. 6), we recorded the magnetization (perturbed (#1)) 
while warming the sample towards the normal state. The 
paramagnetic values for “perturbed (#1)” fall smoothly 
as the temperature approaches Tc, however, this curve 
does not merge with the pristine Mpcc (T) data. Mag¬ 
netization data in the “perturbed {^2)” (open triangles) 
and “perturbed (#3)” (stars) states showed a switching 
towards diamagnetic region after the impulses were ap¬ 
plied at 7K and 6.71 K, respectively. While warming-up, 
the Mpcw curve for “perturbed (#3)” merges with that 
of “perturbed (#2)” at T = 7K, and the two curves to¬ 
gether follow nearly the same path thereafter. It is a mere 
coincidence that the first magnetization data after the 
impulse treatment at 7K in the case of perturbed (#2) 
falls on the warm-up curve of perturbed (#3). Taking 
cue from the observations made in Fig. 3, the magne¬ 
tization after the impulse treatment at 7K could have 
achieved any value ranging from paramagnetic to dia¬ 
magnetic. Some examples of path-dependent Mpcw (T) 
responses in the close proximity of Tc are displayed in 
an inset panel of Fig. 6. In the fourth run, we applied 
the impulse at T = 7.10 K, and recorded the Mpcw {T) 
thereafter while warming up the sample (see Mpcw for 
“perturbed (#4)” in the inset of Fig. 6). This curve can 
be seen to follow a path different from the pristine cool- 



Figure 6: (Color online) Mpcw {T) curves at 77 = 50 Oe 
traced in five different warm-up runs after perturbation by 
an ac field impulse (perturbed (#1) to (#5)). The pristine 
curve (open circles) which was recorded during field-cooling in 
same field is also shown. Magnetization data in the perturbed 
conditions (^^1 to #5) were recorded while warming the sam¬ 
ple to higher T, after perturbing the pristine FC state created 
at T = 7K (perturbed (#1) and perturbed (#2)), T = 6.71 K 
(perturbed (#3)) and T = 7.1 K (inset shows perturbed (#4) 
and (#5)). In all the five cases, the Mpcw (T) curves show 
a variety of paths that can be traversed while warming after 
the perturbation. 



Figure 7: (Color online) Portions of the Mpcc (T) curves 
recorded in 77 = 50 Oe at various amplitudes of vibration 
in a VSM, as indicated. The saturated values of Mpcc at 
T ~ 1.8 K for two of the curves, namely those recorded for 
0.5 mm and 8 mm, can be seen to be different in the inset 
panel. 


down one, as is apparent from the inset panel. When 
we once again applied the impulse at T = 7.10 K i.e., 
the 5*^ attempt, the Mpcw{T) curve (perturbed (#5)) 
nearly traced the path of the pristine Mpcc curve. Thus, 
during the warm-up runs also, one can observe a variety 
and path dependence in M{T) response after creating 
different perturbed states (see “#1” to “#5” in Fig. 6). 


D. Diversity in Mpcc across the region of PME 

We have so far restricted our discussion to the effect 
of a perturbation in the form of an ac field impulse on 
the FC magnetization response (Figs. 3 to 6). It is also 
tempting to ask whether a multiplicity in M (T) similar 
to that observed in Figs. 3 to 6 can also be explored by 
any other form of perturbation, say, via a change in the 
experimental conditions. In this spirit, we show in the 
main panel of Fig. 7, a variety in the Mpcc (T) curves 
obtained at a chosen field of 50 Oe, when the amplitude 
of sample vibration in SVSM was changed from 0.5 mm 
to 8 mm. The following features in Fig. 7 are noteworthy: 

1. The Mpcc {T) curves at different amplitudes tra¬ 
verse different paths, which is very striking. In par¬ 
ticular, the saturated values at the lowest tempera¬ 
ture are significantly different during different runs 
(data shown only for the amplitudes of 0.5 mm and 
8mm in the inset of Fig. 7). Such a difference in 
saturated Mpcc {T) is very similar to the obser¬ 
vations made earlier for the Mpcc{T) responses 
of the pristine and perturbed FC states (cf. inset 
panel of Fig. 3). 

2. The paramagnetic peak can be seen to have high¬ 
est value during the Mpcc (T) run with an am- 
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plitude of 0.5 mm. The peak height was observed 
to steadily decrease as the amplitude was progres¬ 
sively increased from 0.5 mm to 7 mm. 

3. At an amplitude of 8 mm, the Mpcc (T) values can 
be seen to be diamagnetic in the entire temperature 
range 1.8K< T < Tc (cf. the curve for 8mm am¬ 
plitude in the main panel and in the inset panel of 
Fig. 7). 

Figure 7 thus illustrates a rich multiplicity in Mpcc (T) 
getting exposed by the change in the amplitude of sample 
vibration. The similarity between the features emanating 
from Mpcc (T) curves of inset panel in Fig. 7 and those 
of inset panel in Fig. 3 is curious. 

It was observed earlier in Fig. 4(b) that an impulse, 
when applied above about 6.8 K, can switch a given dia¬ 
magnetic value to a paramagnetic value (see perturbed 
(#1) in Fig. 4(b)). Along this line, it will now be in¬ 
teresting to investigate whether such a change can also 
be induced by an impulse in a Mpcc (T) response at 
higher amplitude, say 8 mm, which in pristine conditions 
remained entirely diamagnetic below Tc (in Fig. 7). To 
demonstrate this, we applied an impulse on the pristine 
FC state created (at 8 mm amplitude) at T = 7.05 K in a 
field of 77 = 50 Oe as shown in Fig. 8(a). For comparison, 
we have appended in Fig. 8 the pristine Mpcc (T) curve 
at 0.5 mm amplitude as well. It is interesting to note that 
the diamagnetic value of the pristine Mpcc (T) curve for 
8 mm amplitude changes to a paramagnetic value after 
the impulse treatment at T = 7.05 K. The Mpcc{T) 
curve after perturbation (open triangles) was traced and 
found to be following a different path as shown in the 
main panel of Fig. 8(a). The saturated value of this curve 
(perturbed, 8 mm) at 2K approaches nearer to that of 
pristine Mpcc (T) curve at 0.5 mm as seen in the in¬ 
set of Fig. 8(a). Similar investigations were carried out 
for the pristine Mpcc{T) curve at 7 mm amplitude as 
shown in Fig. 8(b). Here, we performed two different 
runs (independently), after the impact of an impulse at 
T = 7.05 K. In one of the runs (perturbed (#1)), a para¬ 
magnetic magnetization was induced at 7.05 K and the 
Mpcc (T) thereafter follows a certain path (open trian¬ 
gles) , which is similar to that seen for 8 mm amplitude 
in the perturbed state (Fig. 8(a)). During the other run 
(perturbed (#2)), we observed that the magnetization 
data changes from near zero value of the pristine Mpcc 
for 7 mm to a large diamagnetic value after the imposi¬ 
tion of an impulse at 7.05 K, and follows a path shown by 
open squares in Fig. 8(b). We draw an important infer¬ 
ence here that the occurrence of paramagnetic response is 
not specific to a given amplitude of sample vibration, as 
it can be unearthed via a perturbation at all amplitudes 
ranging from 0.5 mm to 8 mm. 

We finally exemplify in the main panel of Fig. 9 an un¬ 
usual feature, viz., a characteristic oscillatory behaviour 
seen in Mpcc (T) response recorded at 0.5 mm ampli¬ 
tude in a higher field, H = 3kOe. This curve exhibits 
some peculiar undulations which override the paramag- 
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Figure 8: (Color online) (a) Mfcc{T) in J7 = 50 Oe, recorded 
at amplitude of 0.5 mm and 8 mm, without imposing a per¬ 
turbation during the cool-down are termed as pristine curves. 
The impulse treatment at 7.05 K switches a diamagnetic value 
of pristine Mpcc (8 mm) into a paramagnetic value, and the 
(perturbed, 8 mm) Mpcc (T) curve was traced on further 
cooling. The inset panel in Fig. 8(a) displays the three curves 
in the range, 2 K < T < 7.5 K. (b) Mpcc (T) at 77 = 50 Oe 
obtained at amplitudes of 0.5 mm and 7 mm, respectively in 
pristine conditions, plotted together with Mpcc (T) obtained 
at 7mm in the perturbed situations {ffl and #2). In per¬ 
turbed (#1) state (open triangles), the Mpcc (T) switches 
from near zero value to a paramagnetic value whereas in 
other case (#2), it switches into the diamagnetic region (open 
squares). 


netic magnetization signal below Tc- It appears as if 
the magnetization shows a tendency to switch between 
the paramagnetic and diamagnetic responses while field¬ 
cooling. Such undulations echo similar oscillatory fea¬ 
tures reported earlier at lower field values (< 14 Oe) 
in another low Tc superconductor, viz., Ca 3 Rh 4 Sni 3 by 
Kulkarni et al. 2^. It may also be added here that 
the oscillatory behaviour in Mpcc (T) is not merely re¬ 
stricted to lower amplitudes and higher fields, but, the 
same were also discernible at higher amplitudes and lower 
fields in 2H-NbSe2. At 7 mm amplitude and in a field of 
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Figure 9: (Color online) A portion of the Mpcc (T) curve 
a,t H — 3kOe recorded at an amplitude of 0.5 mm showing 
some unusual undulations. The inset panel shows a dip (see 
encircled portion) in Mfcc (T) curve at H = 50 Oe, just be¬ 
low the onset of the superconductivity during run with 7 mm 
amplitude. 


H = 50 Oe, the Mpcc (T) response was observed to have 
an oscillatory character in close proximity of Tc as can be 
seen in an expanded portion of the Mpcc (T) in an inset 
panel of Fig. 9. This curve first shows a pronounced dip 
(encircled in the inset of Fig. 9), which can be taken as 
the onset of diamagnetism. Following this dip, there oc¬ 
curs a sharp upturn on further lowering of temperature. 
The Mpcc (T) values eventually cross over to negative 
values on lowering the temperature below T 7.05 K. 


IV. DISCUSSION 

We have come across some amazing metastability ef¬ 
fects in magnetization response on field-cooling a single¬ 
crystal of 2H-NbSe2, viz., occurrence of both paramag¬ 
netic and diamagnetic response in a temperature window 
just below Tc, multi-valuedness in Mpcc {T) at a given 
{II,T) value, oscillatory behaviour in a given Mpcc {T) 
curve, unpredictable switching of Mpcc (T) between dia¬ 
magnetic and paramagnetic values due to an intervention 
by an external perturbation, etc. The mechanism and 
origin behind the said features in 2H-NbSe2 cannot be 
comprehended on the basis of existing theoretical treat¬ 
ments in the literature. The following descriptions form 
the basis of this assertion: 

1. The occurrence of (anomalous) paramagnetic mo¬ 
netization in conventional bulk superconductors [l3| 
had lead to several proposals, amongst which the flux¬ 
trapping model by Koshelev and Larkin received 
significant attention. This model treats the situation 
wherein the magnetic flux gets trapped while field-cooling 
within a (macroscopic or mesoscopic) superconducting 
specimen due to an inhomogenous cooling While 


a (diamagnetic) current flows on the surface to shield 
an external field, associated with the trapped flux in¬ 
side the specimen, there occurs a (paramagnetic) current 
that circulates in the interior of a superconducting sam¬ 
ple . Paramagnetic magnetization response can result 
due to an unbalancing effect of magnetization responses 
arising from the above stated two currents. The inho¬ 
mogeneity in field-cooling, can have a variety, which in 
turn, can lead to a variety in trapped flux, thereby, re¬ 
sulting in multitude of magnetization responses during 
field-cooling experiments. Such a possibility could ratio¬ 
nalize the observation of variety in Mpcc (T) response, 
as reflected in Fig. 7 in the form of vibration amplitude- 
dependent Mpcc (T) data at a given field value. Car¬ 
rying this argument further, the data in Fig. 7 would 
imply that the inhomogeneity in cooling may be lesser 
at higher value of amplitude of vibration (see Mpcc (T) 
curves at 7mm and 8 mm amplitude in Figs. 7 and 8), 
leading to insignificant flux trapping (absence of para¬ 
magnetic signal). However, even in case of Mpcc{T) 
data recorded at higher amplitudes (say 8 mm), that had 
remained entirely diamagnetic below Tc, a diamagnetic 
Mpcc (T) value at a given temperature can be switched 
into a paramagnetic one, when the sample is perturbed 
by an impulse of an ac field (cf. F^. 8(a)). From the 
perspective of flux-trapping model [3l| , this would imply 
that the impulse treatment can drive the system from 
a state of insignificant trapped-flux (diamagnetic) to a 
state of large trapped-flux (paramagnetic), which sounds 
to be unfeasible. Along similar lines, it can be argued 
that the occurrence of a diamagnetic dip just below Tc 
(insignificant trapped-flux) followed by a paramagnetic 
peak (large trapped-flux) in Mpcc (T) (cf. Fig. 9) can¬ 
not be rationalized by this model. Also, the switch¬ 
ing tendency of magnetization between diamagnetic and 
paramagnetic values identified via the characteristics un¬ 
dulations (oscillatory behaviour) in a given Mpcc{T) 
run (Fig. 9) needs comprehension which is beyond the 
scope of flux-trapping model. We further noted two 
different effects of an impulse in two independent runs 
on the Mpcc (T) response obtained at a higher ampli¬ 
tude of 7mm (Fig. 8(b)). In one of the instances (per¬ 
turbed (#1)), the impulse switches a diamagnetic (pris¬ 
tine) Mpcc (T) value into a paramagnetic one while in 
the other (independent) run (#2), it switches the same 
diamagnetic Mpcc (T) value into a larger diamagnetic 
value. Clearly, this suggests that the effects of an impulse 
in the present study are completely random, whereas the 
Koshelev and Larkin’s mathematical treatment [IH to 
a given field and current distribution due to a trapped 
flux ought to be predictable. We therefore rule out that 
this model can explain the metastability observed in 
Mpcc (T) response of 2H-NbSe2. 

2. To understand the thermomagnetic history depen¬ 
dence observed in the magnetization data of a type-II 
superconductor, Clem and Hao [i^ had constructed a 
model which relies on the Bean’s Critical State frame¬ 
work [ 4 ^. Although this model [1^ yields the flux- 
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density profiles that can get set up in a superconductor 
during different protocols usually studied by experimen¬ 
talists, viz., zero field-cooled, field-cooled cool-down and 
field-cooled warm-up modes, it cannot account for para¬ 
magnetic magnetization in the superconducting speci¬ 
men. However, we noted that the effect of an impulse 
applied below T = 6.7 K does not lead to any param¬ 
agnetic magnetization but an enhanced diamagnetism 
(cf. Fig. 5). The multi-valuedeness in (diamagnetic) 
Mpcc{T) data at low temperatures (T < 6.7 K, cf. in¬ 
set panel of Fig. 3, 4(a) and 7 and 8(a) and Fig. 5) 
can be related in an oblique way to a result that em¬ 
anated using the model of Clem and Hao [i^. We re¬ 
call that the flux expulsion by a pinned-superconductor 
when it is field-cooled through its Tc is never complete 
[i^l and therefore the field-cooled magnetization gen¬ 
erally does not conform to an equilibrium magnetiza¬ 
tion. The model by Clem and Hao [i^ indeed shows 
that the Mpcc (T) curve at a constant field depends 
on the pinning strength; a more (less) diamagnetic state 
amounts to less (more) strongly pinned vortex configu¬ 
ration (see Fig. 5 in Ref. [i^). We observe a diversity in 
Mpcc {T) data (see, for example, inset of Fig. 3) that 
closely resembles the pinning-dependent (multi-valued) 
Mpcc {T) curves numerically calculated by Clem and 
Hao (see Fig. 5 in Ref. El)- however, the multiplicity 
seen in present case is in the same 2H-NbSe2 crystal, i.e., 
without changing the amount of quenched-disorder (pin¬ 
ning strength). Comparing the two pictures, viz., the 
Mpcc {T) response of 2H-NbSe2 (an inset of Fig. 3) and 
that shown in Fig. 5 in Ref. |42l| . one may naively sur¬ 
mise that the non-uniqueness in Mpcc {T) in 2H-NbSe2 
could be due to the possibility of coexistence of several 
(metastable) vortex configurations with varying pinning 
strength. The impulse treatment or amplitude-variation 
can lead to the switching of the system from one config¬ 
uration to another, which is reflected by the diversity in 
Mpcc (T) as shown in inset panels of Figs. 3 and 7. In a 
given field-cooled run, the Mpcc {T) values get saturated 
below a certain temperature, Td, at which the applied 
field just equals the lower critical field {Hci{T)). Below 
Tci, i.e., in the Meissner phase, the applied field remains 
smaller than the Hci(T) value and the application of an 
ac field impulse (/loc) superimposed on the (static) ap¬ 
plied field {H ± hac) is hardly expected to produce any 
significant change in Mpcc (T). This is apparent from 
the minimal changes seen in Mpcc (T) below T = 5K 
for the perturbed cases (#2), (#3) and (#4) in Fig. 5. 

3. 2H-NbSe2 system has for long remained a favourite 
compound to explore novel notions in vortex state stud¬ 
ies. In particular, this compound has been very widely in¬ 
vestigated (see Ref. Q and references therein) for order- 
disorder transitions as fingerprinted via anomalous vari¬ 
ations in field/temperature dependencies in critical cur¬ 
rent density, viz., peak effect (PE) phenomenon/second 
magnetization peak (SMP) anomaly in magnetization 
hysteresis loops. Due to supercooling/superheating ef¬ 
fects [i^ that can occur across these transitions, the 


magnetization response is generally found to be history- 
dependent [ij. However, we would like to emphasize 
here that the present results pertaining to metastabil¬ 
ity effects in 2H-NbSe2 have no correlations with the 
metastability and thermomagnetic history-dependence in 
critical current density values seen across the PE/SMP 
transitions. Firstly, most of the magnetization data 
presented here have been recorded at very low fields 
(i.e., H = 50Oe). We have checked that at this field 
value, the present crystal does not display either PE 
in temperature-dependent isofield scans or SMP transi¬ 
tion in isothermal scans at any temperature. Secondly, 
the history-dependent magnetization across the PE/SMP 
has never been reported to generate paramagnetic mag¬ 
netization. 

4. Clearly, all the novel proposals [1, H, S l27l - [30j| that 
have been put forward to understand the origin of para¬ 
magnetic magnetization in high-Tj, superconductors are 
very difficult to visualize in the present context which in¬ 
volves a conventional low-Tc superconductor. Even the 
controlled manipulation of paramagnetic magnetization 
demonstrated recently [2^ in Pb-Co nanocomposites in¬ 
volved the change in the orientation of magnetic moments 
of Co nanoparticles unlike the present situation wherein 
we have observed in a non magnetic superconductor, 2H- 
NbSe 2 , an unpredictable switching of paramagnetic mag¬ 
netization into diamagnetic and vice versa. 

5. We rule out that the present results are a conse¬ 

quence of any specific disorder present in the anisotropic 
2H-NbSe2 crystal as there have been few other reports of 
similar metastability effects (including oscillatory magne¬ 
tization behaviour) in other isotropic low-Tf. supercon- 
ductors, viz., single-crystals of Ca 3 Rh 4 Sni 3 and Nb 
(in the form of a s phe re) [Tst ] . An advantage with a spher¬ 
ical (Nb) crystal is that the demagnetization factor 
remains the same irrespective of the sample orientation 
with respect to field and, hence, yields identical results 
in all orientations. In these studies [l^ however, 
the effects of perturbation on Mpcc {T) response were 
not explored. The novel consequences of perturbation on 
the Mpcc (T) response in the present study have lead 
to exposure of a rich multiplicity in FC magnetization 
and random switching tendency of Mpcc (T) between 
paramagnetic and diamagnetic values, which is indeed an 
advancement over the results reported in Ref. [H, . 

We now focus our attention onto some theoretical 
works [ 32 ], m, [36l - l40j| which attempt to account for the 
occurrence of PME in mesoscopic-sized samples. Based 
on the Ginzburg-Landau theory, the theoretical pictures 
framed exclusively for mesoscopic superconductors fore¬ 
see the occurrence of giant vortex states with multiple 
flux quanta (T > 1) at the onset of surface superconduc¬ 
tivity [^ . Each multi-quantum state has its own field 
domain of existence; at higher field end of each domain, 
the response is diamagnetic, which crosses over to the 
paramagnetic values as the field decreases (see Fig. 2(b) 
in Ref. Figs. 16, 17 and 23 in Ref. E^). It is also pos¬ 
sible that the temperature (or field) sweeping in a given 
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M{T) (or M-H) run is sufficient enough to drive the sys¬ 
tem away from the field/temperature domain of a certain 
configuration (L), which may induce a forced transition 
[ 3 ^ 1^ to the nearest minimum energy configuration 
(with different L value). An examination of the M-H 
plots for different multi-quanta states (L > 1) along with 
the corresponding free energy vs field curves drawn for a 
specific circumstance in Figs. 2 and 3 of Ref. tells us 
that if the transitions occur between L states such that 
the configuration pertaining to a certain L value with the 
lowest free energy prevails, there would be no metasta¬ 
bility effects in the magnetization response and the mag¬ 
netization values as a function of temperature (in isofield 
scans) would be diamagnetic all through. However, if 
metastability is permissible and transitions happen be¬ 
tween different metastable L states, one can encounter 
all sorts of possibilities, i.e., the magnetization values can 
change from a given paramagnetic/diamagnetic value to 
a higher or a lower value, and it can also transit from 
positive to negative values or vice versa. 

To our know ledge, the above mentioned (and other 
related) studies [^, [36l440l| describing the magnetiza¬ 

tion behaviour of multi-quanta states (L > 1) in meso¬ 
scopic samples are the only examples which bear some 
resemblance with our findings. We draw here so me p ar- 
allels between these theoretical findings [H, [H, I36l440l | 
and our experimental results in a sin gle crystal of 2H- 
NbSe 2 , to which the latter results [^, [33, prima 

facie do not apply. As mentioned above, each multi¬ 
quantum state {L > 1) has a certain domain of its ex¬ 
istence wherein the magnetization response can be ei¬ 
ther diamagnetic or paramagnetic depending upon the 
external field [3^. Similarly, in 2H-NbSe2, there exists a 
temperature window (6.7K< T < Tc), where a pertur¬ 
bation applied to the system (at a given {H,T) value) can 
randomly yield the magnetization which is either param¬ 
agnetic or diamagnetic (see Figs. 3, 4(b), 6 and 8). The 
switching of Mpcc {T) due to a perturbation leads to a 
rich diversity in magnetization response in 2H-NbSe2 (cf. 
Figs. 3 to 6 and 8), very similar to characteristics of mag¬ 
netization behaviour anticipated during the transforma¬ 
tion amongst various multi-quanta (L > 1) states [^ . l40l| . 
Additionally, the oscillatory behavior seen in Mpcc (T) 
curve ( Fig . 9) mimics the magnetization response pre¬ 
dicted |33l . |40|| for the transformations amongst multi¬ 
quanta states {L > 1) during the temperature-sweep. In 
an earlier unrelated study, the scanning tunneling mi¬ 
croscopy technique employed by Karapetrov et al. [i^ . 
enabled them to directly observe a coexistence of multi¬ 
quanta vortex states and Abrikosov lattice in an array 
of submicroscopic metal (gold) island embedded in a sin¬ 
gle crystal of 2H-NbSe2. The simultaneous existence of 
multi-quanta state of different vorticities is thus not an 
unrealistic possibility. 

In spite of a strong resemblance of our results with 
the characteristics of multi-quanta states [13, lH, , 

the possibility of occurrence of multi-quanta states in the 
sample investigated here remains a speculation consid¬ 


ering the larger size of this (macroscopic) sample. The 
effect of multi-vorticity in such a macroscopic sample is 
an exotic proposition. Therefore, there is currently a ne¬ 
cessity to look for a more sophisticated model to fully un¬ 
derstand the new findings in 2H-Nb Se9 as well as earlier 
results in CasRh^Sni^ fEA l and Nb flaJ . 

In the end, we rule out the possibility of any exper¬ 
imental artefact arising due to change in amplitude of 
vibration considering the following important points: (i) 
The factor which could be affected most by changing 
the sample vibration is the magnetic field-inhomogeneity 
which, in the present case, is found to be negligible ( 
~ 10“^ Oe). (ii) We did not register any change in the 
magnetization on varying the amplitude at a given {H,T) 
value in the case of a standard Indium sample, (iii) We 
have also shown that the PME like attribute in 2H-NbSe2 
is not limited to only smaller amplitudes but the same 
can also be induced at higher amplitudes (8 mm), after 
the impulse treatment. 


V. CONCLUSION 

We have investigated in detail a single crystal speci¬ 
men of 2H-NbSe2 via dc magnetization measurements. 
The new revelations in this compound correspond to the 
observation of an anomalous paramagnetic (PME like) 
signal(s) below Tc and some peculiar consequences of ap¬ 
plying a perturbation at various temperatures on a given 
field-cooled state. Based on the effects of perturbation, 
two distinct temperature intervals have been identified. 
Across a certain temperature range lying just below Tc, 
the PME - like signal can be manipulated by applying 
an external perturbation, which unpredictably switches 
the Mpcc (T) response from a given para/diamagnetic 
value into a different para/diamagnetic value. On the 
other hand, when the system is perturbed away from 
Tc (below about T = 6.7K), the effect of perturbation 
is of one kind, i.e., it only enhances the diamagnetism 
such that no paramagnetic signal is seen. The saturated 
(diamagnetic) Mpcc (T) value at the lowest temperature 
is found to be influenced by the effect of perturbation 
applied closer to the region of PME. The said features 
in Mpcc (T) data in a single crystal of 2H-NbSe2 co¬ 
incidentally bear similarities with the magnetization re¬ 
sponse predicted [13, [H, [36[ - [43| for mesoscopic samples, 
wherein, a possibility of nucleation of multi-quanta states 
{L > 1) have been discussed in the literature. We believe 
the present results in a bulk 2H-NbSe2 crystal cannot be 
explained by any of the existing mathematical treatments 
based on prescriptions of Bean’s Critical State model [131 
reported in the literature and call for a new theoreti¬ 
cal framework. Though, we have presented experimen¬ 
tal results pertaining to a given orientation of magnetic 
field vis. a vis. crystallographic direction {H\\ab) of 
hexagonal 2H-NbSe2, we have also carried out magneti¬ 
zation measurements in more often studied orientation, 
viz., H\\c, in 2H-NbSe2 and obtained nearly identical 
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features (data not included in the present report). We 
may therefore state that assertions made in this report 
are not specific to the anisotropic nature of the intrinsic 
superconducting parameters of 2H-NbSe2. 
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